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The effect of lielical Dirac states on surface phonons in a topological insulators is investigated. 
Their coupling is derived in the continuum limit by assuming displacement dependent Dirac cones. 
The resulting renormalisation of sound velocity and attenuation and its dependence on chemical 
potential and wave vector is calculated. At finite wave vectors a Kohn anomaly in the renormalized 
phonon frequency is caused by intraband-transitions. It appears at wave vectors q < 2kF due 
to a lack of backscattering for helical Dirac electrons. The wave vector and chemical potential 
dependence of this anomaly is calculated. 
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I. INTRODUCTION 

Strong topological insulators^ have metallic surface 
electron states which are protected by time reversal 
invariance^. One of the most prominent examples is 
Bi2Se3 which has a single Dirac cone with a linear 
dispersion at the zone center of the surface Brillouin 
zone. Due to the strong spin orbit coupling the spin 
is locked perpendicular to the wave vector for each 
surface state. These features have been confirmed by 
photoemission^^. The chemical potential may be tuned 
by surface or bulk doping* or application of a gate 
voltage^. Apart from photoemission it has been diffi- 
cult to find the signature of surface states in transport 
results which may be masked by bulk contributions^. 
It is therefore desirable to use methods that probe 
naturally only the surface states. One possible method 
is to study the propagation of surface acoustic waves 
(SAW) by ultrasonic means. This method has been 
used successfully for the conventional two dimensional 
electron gas in GaAs heterostructures with parabolic 
dispersion^'^- and it has also been proposed for studying 
magnetoconductance of Dirac cone quasiparticles in 
grapheneii. 

Here we investigate the surface acoustic phonon 
propagation in topological insulators, using the con- 
tinuum limit for the surface conduction state and the 
phonon modes. This means a linear dispersion in both 
cases given by Fermi velocity vp and sound velocity Vs 
respectively. Experimentally so far only bulk optical 
phonons have been investigated in Bi2Se3i^ii^. The 
renormalisation of surface acoustic phonons by helical 
Dirac electrons should carry a signature of the latter, 
i.e., their linear density of states around the Dirac point, 
the linear dependence of Fermi wave number on chemical 
potential and the suppressed backscattering due the spin 
locking perpendicular to the wave vector. The latter 
should play an important role at finite phonon wave 
vector where the phonon frequency is renormalised by 
both intraband and interband contributions. 



In Sect. |TI] we give the basic properties of the he- 
lical Dirac fcrmions. In Sect. IIIII their coupling to 
surface phonons is derived and in Sect. IIVI the renor- 
malised sound velocity, attenuation and Kohn anomalies 
in phonon frequency are calculated. Finally Sect. IVl gives 
the conclusions. 



II. HELICAL DIRAC FERMION MODEL 

The protected surface states of a 3D strong topologi- 
cal insulator form a single Dirac cone located at a time 
reversal invariant point in the Brillouin zone. Their low 
energy (as compared to the bulk valence band gap) effec- 
tive Dirac Hamiltonian may be written aa^^^— 



Hq = ivp I (frtl;\r)(T ■ (z x V)V'(r) 



(1) 



Here vp the Fermi velocity, (t,z are the electron spin 
and the fixed surface normal respectively. The quantisa- 
tion axis for the electron spin is also chosen parallel to 
the surface normal. The free massless Dirac electrons are 
represented by the two component spinors 



V^(r) 



1 



d re 



Ckt 
Cki 



(2) 



Adding a term including the chemical potential /i the 
diagonalsation leads to the Dirac cone energies 



Ckr = ±WF|k| - /i 



(3) 



where |k| = {k'^ + ky)^ and t — ±1 is the helical pseu- 
dospin index denoting the postive or negative spin helic- 
ity of the eigenstates. They are obtained from the free 
electron states by the unitary transformation matrix 



W : 



le 



(4) 



Which leads to the creation and annihilation operators 
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FIG. 1. Schematic view of adiabatic modification of helical 
surface states by a surface phonon in the xz- plane. The 
surface normal n(j:) = z + Sn{x) giving the local orientation 
of the Dirac cone is periodically modulated by the surface 
displacement amplitude u(x). 



for helical eigenstates 7^^ according to 



T le' - 4, ) 



kj,/ 



7k± = -^(e' 2 ckt ± ?e * 2 ck^) 



(5) 



Here the rotation angle to helical eigenstates is given 
by tan ??k — k^/ky. The existence of surface helical states 
in the Dirac cone has been proven most clearly for the 
Bi2Sc3 compound by ARPES experiments^. As men- 
tioned in the introduction further evidence for their phys- 
ical effects should be primarily sought with other surface 
specific probes. 



III. COUPLING OF SURFACE PHONON 
MODES TO THE DIRAC CONE 



One possibility is the investigation of long wave 
length surface phonons, in particular surface acoustic 
waves (SAW). The latter has been very successfully 
used for the study of the 2D electron gas with parabolic 
dispersion found in GaAs heterostructures^'^^ and it 
has also been proposed as a method to investigate the 
chiral Dirac electrons in graphene monolayer; 
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particular the effect of Landau quantisation. Similar 
considerations may be fruitful for the surface Dirac 
states of 3D topological insulators. In the previous 
cases a phenomcnological approach was used where the 
modification of surface phonon properties is expressed 
through the conductivity tensor. Here we prefer a more 
microscopic approach to investigate the effect of the 
unconventional helical spin polarisation of surface states 
in the topological insulator. 

We first derive a suitable low energy and long wave 
length effective Hamiltonian for the coupling of surface 
phonons and helical Dirac states. This will be done by us- 
ing an adiabatic generalisation of the unperturbed Hamil- 
tonian in Eq. ([T]) valid for <C wf where is the sound 



velocity. This can be achieved by replacing the global 
surface normal z by the local normal n(r). Its depen- 
dence on surface position r = {x,y) is determined by the 
amplitude of the surface phonon mode. This leads to the 
modified Hamiltonian 



H = ivp 



j dPvi:\r)(T ■ (n(r) x V)V'(r) 



(6) 



For small phonon wave vectors (g ^ fco = 7r/a, a 
~ lattice constant) and amplitudes we may expand 
n(r) = z + (5n(r) where (5n(r) is the local deviation from 
the global surface normal due to the phonon mode. The 
geometry of this deformation is schematically shown in 
Fig. [TJ For simplicity we consider only displacements 
perpendicular to the surface. In the long wave length 
limit surface acoustic phonons also have a parallel com- 
ponent enforced by the stress free surface condition, but 
this may be taken into account in the end by a simple 
numerical factor as discussed in Appendix \^ Figure [1] 
illustrates the origin of the Dirac electron-phonon 
coupling. The continuous long wave length deformation 
of the surface leads to a periodic tilting of the surface 
normal. The local quasiparticle states have Dirac cones 
alligned with the local surface normal. We consider only 
the case of one cone in the surface Brillouin zone. In the 
continuum limit without inter-cone scattering the effect 
of several cones would simply be additive. Expanding 
with respect to the tilt angle one obtains a coupling 
of the unperturbed quasiparticle states to the surface 
phonon coordinates. 

In the following we derive the relevant coupling Hamil- 
tonian. Denoting the displacement vector by u(r) the 
local direction of the surface normal is given by 



n(r) 



1 



(1 



(7) 



Approximation to first order in the displacement ampli- 
tude leads to 



fi(r) ~ z + (5n(r) 
(5fi(r) = — u'^x 



(8) 



The first term is the global surface normal z and the sec- 
ond the local correction due to the surface displacement 
where u'^ — duz{x)/dx is the displacement gradient par- 
allel to the phonon wave vector q = (?q where q = x. 
Expanding in phonon coordinates we get 



(9) 



Inserting n(r) into the Hamiltonian of Eq. © we can split 
into unperturbed and Dirac electron-phonon Hf,p cou- 
pling term. Rearrangement of terms in the latter finally 
leads to 



(10) 



Hep = -^X!fq[(k X q) -Z] 
kq 



^i+qa,Vk](6q + 6Lq) 
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with the couphng constant defined by 



5q = (1[ 



2Mujc 



-(^) (11) 



or g„ = 



Here Wq — VgQ is the frequency and Vs the velocity of 
long wave length surface phonons. It is related to the 
elastic constant Cg — pvl with p = M/Vc denoting the 
density and Af, Vc the atomic mass and volume per 
atom respectively. Because [(k x q) • z] = sin6'kq where 
^kq = Z(k. q) the electron- phonon coupling in Eq. PT|) 
vanishes for forward scattering (0kq = 0) and backward 
scattering (0kq = ±7r). In the following we simplify the 
notation to 0k = ^kq because q = ic is fixed. 

To perform useful calculations the coupling term has 
to be transformed to the helical Dirac eigenstates of Hq 
given in Eq. ([5]). This is accomplished by applying the 
unitary transfromation W to Hf,p. We then obtain 

1 



"^^]^Z!fq[(kx q) 

kq 



zli^UqfkqVfkKbq 



7k+ 
7k- 



Tkq — rkqTo + rkq^a; 



+ 6Lq) 



(12) 



Here tq, Tx are unit and Pauli matrix respectively in the 
space of helical eigenstates. The intra- (F) and inter- 
(r) band transition matrix elements between opposite 
and equal helicity states respectively are given by {k± = 
kx i iky etc.) 

■{k + q) + k 



Fkq = k sin ^k 



Tkq = k sin 6k 



|k + q|fc 

(fc + g)+fc_ 



|k + q|fc 



1 



1 



(13) 



The momentum dependent part in the parenthesis may 
be written as 

{k + q)+k_ fc + k • q + i{k x q) • z 



|k + q|A: 



(fc2 



ge 



2k- 



(fc2 



2kqcos ^^k)' 



(14) 



In the long wave length case g/fc <C 1 we can approximate 
Tkq = iq sin^ Ok 

fkq = 2fcsin6'fc (15) 

This limit corresponds to the physically relevant case for 
ultrasonic SAW provided that q ^ kp — p/vp which 
holds except for the chemical potential very close to the 
Dirac point. 



IV. SURFACE PHONON RENORMALISATION 
BY DIRAC PARTICLE-HOLE EXCITATIONS 

In this section we calculate the phonon self energy and 
derive the renormalized sound velocity and phonon fre- 
quency and in particular the damping due to excitation 



of surface state conduction electrons. The dependence on 
wave vector and chemical potential of these quantities is 
in principle accessible by experiment. 



A. Phonon self energies 

The frequency ujqs of surface phonons renormalised by 
coupling to helical Dirac electrons is obtained from the 
Dyson equation as 



-qs • "-^qsns(q, Wqs) 

ns(q,Wqs) = n(q,Wqs) -hn(q,(I)q5) 



(16) 



where ns(q, cjqs) is the total self energy due to surface 
electrons and H, H denote the intra- and inter-band con- 
tributions respectively. The effects of the background 
bulk excitations across the valence-conduction band gap 
are assumed to be already included in the unrenormalised 
phonon frequency Wqg . In second order perturbation the- 
ory with respect to the coupling constant gq the contri- 
bution of surface states to the self energy is given by 
(a* > 0): 



2„2Y^iT.. |2^( /^ - ^F|k + q|) - 0{p. - vf\^\) 
a; + t;F(|k + q| - M) + irj 



n(q,^) = 44^|rkq 



n(q,c.) = 445]|fkq| 



-e{vF\\^ + (i\- p) 

w + WF(|k + q| + |k|) + iry 



(17) 



The excitation processes contributing to the self energies 
are indicated in the inset of Fig. [2] We do not evalu- 
ate these expressions in the general case but restrict to 
a; — > appropriate in the adiabatic limit <C fF- Defin- 
ing n(q, w) = u|n(7qn(q, oj) etc. we can approximate for 
phonons with q/kp <C 1 according to 

ii{ci) = -Y^\Tkq\^5{^,-vp\\.\) 



ft(q) = -E 



kql 



2?;F|k 



■e{vp\\<i\- p.) 



(18) 



Using the long wave length expressions for the matrix 
elements given in Eq. (jl5p we evaluate these expressions 
as 



n(g) = --D{^,)q^ 

m = -\{^)D{p){ki-ki) 



(19) 



where k^. = p-c/vp is a cutoff which corresponds ap- 
proximately to the wave vector where the surface states 
merge into the bulk bands. Here /ic is one half of 
the conduction-valence band gap. Furthermore D{p) — 
is the density of states at the chemical potential 

with A — a? denoting the surface unit cell area. Obvi- 
ously for q/kp <§; 1 the intraband contribution is negli- 
gible compared to the interband part because the former 
involves only few excitation processes around p. 
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FIG. 2. Renormalised phonon frequency change As(q) — 
{Ldqs — ajqs)/tjqs for /i//ic ~ kp/kc = 0.5. Interband contri- 
bution is monotonous while intraband part exhibits a Kohn 
anomaly for q/kc — q/2kF < 1 which is also visible in the 
change of total renormalised phonon frequency. The inset 
shows a schematic view of intra- and interband excitation 
processes. 



B. Sound velocity correction and attenuation of 
surface acoustic waves 

Then, using Eqs. (|16ll9p one finally arrives at the 
surface phonon velocity correction (p > 0) defined by 
Avs — {vs — Vs)/vs- We obtain 

where fco = ti'/o defines the dimension of the surface 
Brillouin zone. We can also vifrite Afw^ = CsVc because 
Vs = (p/cs)"^ ■ Equation (|20l) is a constant velocity shift, 
independent of q, but it depends on the Fermi level fi. 
It may therefore be identified as surface contribution by 
changing the chemical potential e.g. by application 
of a gate voltage. In fact it vanishes when the Fermi 
level reaches the conduction band aX ji = pc because no 
possible surface interband excitation processes are left 
in this case. 

A further telling effect of surface Dirac electrons may 
be found in the intrinsic attenuation coefficient due ab- 
sorption of surface phonons by particle-hole excitations. 
For that purpose the imaginary part of the self energy at 
small but finite frequencies ui = VsQ has to be evaluated. 
It is obvious from Eq. (IT71) and the inset in Fig. [5] that 
only intraband processes can contribute to the phonon 
absorption because <C /i. They are given by the 
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FIG. 3. The evolution of Kohn anomaly in the phonon fre- 
quency change As(q) = (ojqs — i^qs)/tt)qs with chemical poten- 
tial or Fermi wave vector /i//ic = kp/kc- For increasing ^{kp) 
the intraband contribution and associated anomaly grows and 
shifts to larger momentum. The values for q — >■ show the 
dependence of sound velocity change Aiig/Aj on the chemical 
potential. 

imaginary part [x = q/k) 

n"(q,a;) = iFkql cos 61k 

k 

■5{li — VFk)5{ijj + XV pk COS 6\s) (21) 
For oj — VsQ this may be evaluated leading to 

n"(g)=A4%2[l-(^)V]^ (22) 

Vp VF VF 

The attenuation rate per length is defined as a^s = 
jqs/fiVs where jqs — Vpgq^" {(i) and in leading order 
of q is finally given by {A — a'^) 

This is an intrinsic damping given by phonon absorption 
due to particle-hole pair creation which does not depend 
on the mean free path of surface electrons. Its charac- 
teristic features are the dependence on wave number 
(since to = Vgq) and the linear dependence on the chem- 
ical potential or surface density of states. The velocity 
shift and attenuation of surface phonons in topological 
insulators can possibly be observed in ultrasonic SAW 
experiments if the change of /i around the Dirac point 
can be achieved by a gate voltage tuning. 

C. Surface phonon Kohn anomaly for helical Dirac 
electrons 

In this section we consider the renormalisation of sur- 
face phonon frequency at finite wave vector compara- 
ble to the Fermi wave vector hp = fJ-fvp- Because the 
self energy will depend on the wave vector one expects 
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a q-dependent renormalisation of the phonon frequency 
which wUl be noticable in particular around |q| ~ 2kp. 
In analogy to the case of bulk bands this may be called 
a Kohn anomaly caused by helical Dirac electrons. Since 
kp <^ ko — n/a the background (bare) phonon disper- 
sion ujqs is still linear in |q| in this regime. We consider 
the normalised change in the phonon frequency given by 
As(q) = {ujqs — Wqs)/a;qi. It may be calculated from the 
static part {uj 0) of the self energies since Vg/vp ^ 1- 
We assume q = gx and define dimensionless quantities 
p, — n/nc, q = q/fcc, k = k/fcc- After some algebra we 
then obtain (/i > 0) 

A,(q) = A°(0(q) + 0(q)) 

0(q) = - / dk deP sin^ 9 
Jo Jo 

k + qcos9. 9{\k - fl\) - 9{\k + q| - A) 



[1 



|k + q| 

r 1 /'27r 

[q) ^ - I dk d9k^ aiv? i 
'o Jo 



|k + q|-|k| 



A: + (7Cos6' 6'(|k + q| - /i) 
|k + q| |k + q| + |k| 



(24) 



where |k + q| — {k^ + q + 2kqcos9)^ and 9 = 9^^ 
is the angle between the wave vectors. The momentum 
integration extends up to the cutoff wave number k — kc 
or k — 1 where the helical surface states merge with 
the bulk bands. The scale of the q-dependent phonon 
frequency renormalisation in dimensionless units is given 

by 



A 



(25) 



It also determines the scale of SAW velocity change 
in Eq. ^ with Av^ = (7r/3)A°. We note that the 
integrands in Eq. (|24p vanish for 9 — t: where the phonon 
couples to back-scattering processes k^k + q= — k. 
This arises because helical states with k and -k have 
opposite spin and the phonon cannot couple to spin-flip 
processes according to the Hamiltonian in Eq. ((TT|) . In 
the usual bulk conduction bands where the spins are not 
coupled to the momentum system it is the backscattering 
intraband processes that lead to the 'Kohn anomaly' of 
renormalised phonons at a wave vector q = 2kp. It may 
be more or less pronounced depending on dimensionality 
and nesting properties of the bandsii. 

In the present case we have both interband and 
intraband contributions to the renormalised phonon 
frequency. The resulting change in the velocity is shown 
in Fig. [5] with both contributions shown in addition 
to the total change for a chemical potential fi — 0.5 
(kp = 0.5fcc), i-e. halfway in the upper Dirac cone. 
The absolute value of the interband contribution shows 
a monotonously decreasing behaviour as function of 



q. On the other hand the magnitude of the intraband 
contribution has a pronounced maximum. Because of 
the suppression of backscattering it is, however, at an 
intermediate wave vector q < 2kp = kc between forward 
(q=0) and backscattering {q — 2kp). As a result the 
total renormalised velocity change shows a resulting 
Kohn anomaly appearing a,t q < 2kp. 

This anomaly in the surface phonon frequency change 
depends on the chemical potential or Fermi wave vec- 
tor. It is shown in Fig. [2] for three different cases. When 
/i = there is no intra-band contribution and one has 
monotonous q-dependence. For increasing /i the Kohn 
anomaly progressively appears and shifts to higher q- 
values. However the maximum effect is always seen at 
q < 2kp due to the lack of backscattering. The absolute 
magnitude of the anomaly is determined by Eq. (P5|) . Us- 
ing parameter estimates for Bi2Se3 according to /ic — 0.3 
eV, CgVc — 4.68 eV and /cc/fco — 0.1 we get the dimen- 
sionless anomaly scale A^ ~ 1.6- 10~^. It is largely deter- 
mined by the small value of (/cc/fco)^ — 0.01 which is the 
small relative area of the Dirac cone at /i = /ic compared 
to the surface Brillouin zonei^. Although A^ is not very 
large it is perfectly sufficient for the observation of elastic 
constant anomalies (as function of n) since for the latter 
relative changes can be measured with a precision up to 
10~^. On the other hand observing a Kohn anomaly at 
this scale is a challenge. In the continuum Dirac model 
there is no nesting involved in the 2D Fermi surface. 
Therefore there is no singular enhancement of the Kohn 
anomaly. However if one includes warping effectsi^ii^ 
due to higher order terms in the surface cone dispersion 
which are important at larger chemical potential nest- 
ing effects appear between various symmetry k- points 
(not connected by backscattering). This could lead to a 
considerable enhancement of the Kohn anomaly at these 
nesting vectors. 



V. CONCLUSION 

In this work we investigated the effect of helical Dirac 
electrons in topological insulators on the propagation of 
surface phonons. Their coupling to the quasiparticle-hole 
excitations has been derived within a continuum model 
assuming a displacement dependent dispersion. The 
coupling to Dirac cone quasiparticles renormalizes the 
sound velocity and leads to intrinsic damping of surface 
acoustic waves. The size of these renormalisation effects 
depend directly on the chemical potential and may 
therefore be checked e.g. by application of a gate voltage 
or opening a gap through magnetic doping^ in surface 
ultrasonic experiments. Furthermore at finite wave 
vectors the surface phonon dispersion exhibits Kohn 
anomalies. Because of the vanishing backscattering due 
to helical spin polarization they appear at wave vectors 
q < 2kp and their magnitude depends on the position 
of the chemical potential. Due to the absence of nesting 
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effect in the continuum limit they show no singular 
q-dependence. The latter may change if warping effects 
are included for higher chemical potential. It will be 
an experimental challenge to identify these predicted 
anomalies which may also be tuned by application of 
gate voltage. 
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Appendix A 

Our discussions of surface acoustic phonon coupling 
were simplified by ignoring the longitudinal displacement 
field that exists for surface waves due to the stress free 



boundary conditions in the limit q — )• 0. We briefly dis- 
cuss here the modifications caused by the inclusion of the 
longitudinal displacement part {x) . In its presence the 
modification to the surface normal in Eq. (|5]) is changed 
to 

v' u" 

Sh{r) ^ + 2^^]x = -<[1 - 27.]x (Al) 

where prime and double prime denote first and second 
derivatives with respect to x (propagation direction) re- 
spectively. Evaluation of the derivatives for a Rayleigh- 
type long wave length surface phonon leads to 



IS — 2 I ( 

(2-a-2(i-a'(i-^a^ 

Here vi^t are longitudinal and transverse bulk sound ve- 
locities and ^ = Vs/vt where Vs is the surface phonon 
sound velocity. Thus inclusion of the longitudinal dis- 
placement field merely modifies the electron phonon cou- 
pling constant in Eq. (fTTj) by a factor (1 — 27^). 
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